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Abstract 

We analyse the entanglement generation in a one dimensional scattering process. The two col- 
liding particles have a Gaussian wave function and interact by hard-core repulsion. In our analysis 
results on the entanglement of two mode Gaussian states are used. The produced entanglement 
depends in a non-obvious way on the parameters ratio of masses and initial widths. The asymp- 
totic wave function of the two particles and its associated ellipse yield additional geometric insight 
into these conditions. The difference to the quantitative analysis of the amount of entanglement 
generated by beam splitters with squeezed light is discussed. 
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I. INTRODUCTION 



In the last decade the study of entanglement has become one of the major topics in 
the flourishing field of quantum information theory (for introductions into this topic see 
. Among other things, it has become clear that entanglement is an essential re- 
source for many desirable operations in quantum information theory, for example quantum 
teleportation |3j] or some protocols in quantum cryptography [4j. For a detailed understand- 
ing of entanglement it is necessary to quantify it. For pure states G this is fairly 
straightforward. We consider the reduced density operator in 7i\ (tracing out system 2) 

p 1 := Tr 2 (|^> <^|) (1) 

and define the entropy of entanglement E {\ip)) of the pure state \ip) as the von Neumann 
entropy of this reduced density operator 

E := S (p 1 ) with S (p) := -Tr (p log 2 p) . (2) 

The quantity E has the nice property that it (asymptotically) describes the number 
of maximally entangled singlet (or Bell) states of 2 qubits that can be obtained from many 
copies of the quantum state when it is shared by two parties and each is allowed to 
operate on his part only jsj]. The quantification of entanglement for mixed states is a much 
more difficult business which however does not concern us in this paper, as we are only 
dealing with pure states. 

In this paper we study entanglement generation in quantum systems with a continuous 
degree of freedom, in particular a scattering process which generates correlations in the 
two-particle Schrodinger wave function. Continuous variable quantum systems are also an 
interesting topic in quantum information theory (cf. the recent review |6|). For instance, 
continuous versions of quantum teleportation p] and quantum cryptography [81 can be 
found. 

In particular we will be dealing with Gaussian states, since properties of these states 
can be often obtained in an analytic fashion such that there is also an extensive theoretical 
.itetatute about these states Q Q. Gaussian wave functions ate also a qui te natuta, 
assumption. We will investigate how the amount of generated entanglement depends on the 
ratio of the widths of the incoming particle wave functions. 
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In quantum optics the entanglement production for photons can be realized experimen- 



tal 
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and has been extensively studied, e. g. for squeezed states hitting a beam splitter 



121 ] . In contrast to that, there are only a few studies of entanglement production for 
scattering of two particles |l3t 1 141 . This is surprising as scattering theory itself is a major 
topic in quantum mechanics. In [1^ entanglement production for scattering of two Gaussian 
particles with hard core repulsion (cf. the potential in eq. (j3J)) was studied, and the present 
paper can be understood as a generalization of this paper which yields additional physical 
understanding. For, in only the special case was treated, when the two particles have 
equal masses and equal initial widths of their wave functions. It was then found that only 
transient entanglement, i. e. non-vanishing entanglement during collision, can be produced. 
For a more general repulsive potential only a rather small quantity of asymptotic (perma- 
nent) entanglement could be obtained with the initial conditions of one can 
find calculations of entanglement in spin gases, where the position coordinates are treated 
classically. 

Our calculation looks quite similar to the calculation of the entanglement that is generated 
in a beam splitter when at least one input state is squeezed. However, the decisive difference 
is that the action of the beam splitter is a rotation that mixes the quadrature amplitudes 
of different modes, whereas the scattering process is a reflection with respect to the non- 
orthogonal coordinate system given by the relative coordinates. 

We now outline the organization of this paper. In section|n]we state the scattering process 
and review the exact solution of the time-dependent Schrodinger equation. In section IIHI 
we calculate the asymptotic entanglement of the two particles after the scattering. Here 
we use extensively the covariance matrix of Gaussian states. A geometric interpretation of 
the results about entanglement is given in section IIVI where we study the wave function of 
the two particles after scattering. Therewith we can also understand in detail under which 
conditions significant entanglement is produced - and also when no entanglement at all is 
generated. In section El we sum up the obtained results and point out some possibilities for 
further research. 
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II. THE SCATTERING PROCESS AND ITS EXACT SOLUTION 



Here we consider the quantum mechanical scattering process of two particles in one 
spatial dimension which interact via the potential V(x\ — a^). As the scattering potential 
shall model hard-core repulsion, it is given by 







x 2 \ 


> a 


( oo , 


\x\ — 


X 2 \ 


< a 



V (\ Xl - x 2 \) := { I -| (3) 

In contrast to we consider general mass ratios of the two particles. Hence the Hamilto- 
nian for the corresponding one dimensional scattering process is 

H =A + ^k +v( ^- xA)= A + A +vM • (4) 

where we have defined the momenta 

p s := pi +p 2 , 

Pr ■= Hi Pi ~ Hi P2 (5) 
and center of mass and relative coordinate x s and x r , respectively, 

x s = Hi x i+ H2 x 2 , 

x r = xi - x 2 , (6) 

using the mass fractions 

jtti :=mi/(mi + m 2 ) , fi 2 ■= ~m 2 / (rnx + m 2 ) . (7) 

Besides, in eq. (jlj) we use the masses M s := (m 1 + m 2 )/2 and M r := (mi m 2 )/(m 1 + m 2 ). 
We thus see that the Hamiltonian (jlj) decouples in the coordinates (x s , x r ). 

Our starting configuration at t = are two Gaussian wave packets with widths af and 
of. Let the mean value of particle 1 be located at Q\ ^> a, that of particle 2 at — Q 2 —a. 
The momentum of the first particle 1 shall be —K and that of the second particle K. We 
first define a state |/o) with wave function 

k(x u x 2 ) = <p G {x x ; Qi, —K, of) <p G (^2 ; -Q2, K, of) , (8) 

where the Gaussian (located at Q and with momentum K) is defined as 

0G ( ; Q, K, o 2 ) := a(o 2 ) exp (iKx) exp y 2 ^ 2 J , a{o 2 ) := —^=—^ . (9) 



Note that |/o) itself does not define a physical initial condition for any finite Q, since the 
Gaussians always overlap. In a more formal approach the two wave packets have to be 
starting from an "infinite distance" jl9l]. 

Now we would like to find a solution \ip t ) E 7i of the time-dependent Schrodinger equation 
with the Hamiltonian H in eq. (@J) that coincides for t — > — oo approximatively with \f t ) : = 
exp(-iH Q t)\fo), where 

h ^A + A ■ (io) 

is the free Hamiltonian. In the relative coordinate x r we have to satisfy the following 
(boundary) conditions: 

(A) For all x T E [—a, a] we have ipt(x s , x T ) = , Wt EM.. 

(B) The solution ipt(x s , x r ) is continuous in x r at the boundary of the potential x r = a. 

(C) For x r G" [—a, a] the solution i/jt(%r, x s ) obeys the free Schrodinger equation with Hamil- 
tonian Hq. 

We now construct the solution as follows. For an arbitrary wave function <p(x r ) we 
introduce the unitary operator P a that reflects the relative coordinate at 

(P a \<i>))(x r ):=<j>(2a-x r ) . (11) 

It can be easily seen that the operator 1 ® P a commutes with the free Hamiltonian (|l()jl 

[1 ® P a , £T ] = . (12) 

Here it is important to keep in mind that the tensor product structure refers to relative 
coordinates and not to the particle Hilbert spaces. Let us define the state 

|<7t):=(l®P«)|/ t > • (13) 

Because of the commutation relation (fT2*|) the evolution t \— > \g t ) is also a solution of the 
free Schrodinger equation. Now the difference \ft) — \gt) satisfies condition (B). In order to 
enforce condition (A) we use a Heaviside function, and thus the general solution that we 
were looking for can be written as 

ipt{x s , x r ) = ( f t (x s , x r ) - g t (x s , x r ) ) 6(x r - a) . (14) 
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Since we have 

lim |||tk)-|/t>||=0 (15) 

t— »— oo 

and 

lim || |^)- =0 (16) 

we consider the mapping \f t ) \— > \g t ) as the scattering process. 

It is important to note that due to the hard-core repulsion particle 1 always stays to 
the right of particles 2. In other words, there is only reflection and no transmission for this 
particular potential. 



III. ENTANGLEMENT BETWEEN THE TWO PARTICLES AFTER SCATTER- 
ING 

In order to calculate the entanglement between the particles after the scattering, we have 
to consider solution (pPfj) for times t —>■ oo. Since eq. (jlfij) holds and the free evolution 
exp(—iH t) does not change the entanglement, we have 

pmE(\iH)) = E(\go)), (17) 
where \g Q ) = (1 Cg> R) |/ ) (cf. eq. ()13|)). It is very important to note that the state \g ) is a 



Gaussian state 
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2Cj For, the initial state |/o) is obviously Gaussian and the reflection in 
relative coordinates does not change the Gaussian nature of the state. In order to obtain 
the entanglement of \go), we associate the following reduced density operator to the state 
\9o) 

p* Ed :=Tr 2 (\g }(g \) (18) 

The reduced density operator of a Gaussian state is also Gaussian. 

The entanglement of a Gaussian state can be determined from its covariance matrix that 
we now introduce. The four standard canonical operators of our two particle system are 

R := (£ ls pi, x 2 , p 2 ) T ■ (19) 

The 4 by 4 covariance matrix of a density operator p is then given as 

a l3 = 1 Tr (p {Rt, Rj}) - Tr (p Rfi Tr (p R,) , i, j G {1, . . . , 4} , (20) 
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where { , } denotes the anticommutator. We recall the following fact (cf. e. g. P|): a linear 
change of canonical operators with a 4 by 4 symplectic matrix F and a displacement D , 
i. e. 

R=FR+D, (21) 
implies for the covariance matrix a of the new canonical operators R 

a = F a F T . 



(22) 



In order to obtain the covariance matrix of the state 



p' :=\g )(go\ = U prf with p := |/ ) </ | , U := 1 ® P a , (23) 

we start with the covariance matrix for p = \f ) (f \ and the standard operator coordinates 
R of eq. (HH 

/ of/2 \ 



a := Tr (RR T p) - Tr (Rp) Tr (R T p) 





1/(2 of) 
of/2 



(24) 



V 1/(2 of)/ 

as can be easily calculated from formulas (jSJ) and Q- The covariance matrix for p' in eq. (|23|) 
in standard canonical coordinates R can be transformed as (using eq. 



a[ 3 = Tr (p' IU R 3 ) - Tr (p' R t ) Tr (p' Rj) 

= Tr (p (f/ t Ri U) (f/ f ^ C/)) - Tr (p (U ] R { U)) Tr (p (f/ t ^ U)) . (25) 

We now seek a symplectic transformation which relates the new operators U* R U to the 
standard operators R. Then we can use formula (|22j) and obtain a 1 from a in (|24j) . Since 
the operator U — 1 (g) P a acts on center of mass and relative coordinate, we define the vector 
S of canonical operators in this frame as 



S (^-s; Ps; ^ri Pr) 

As can be easily seen from eqs. © and ©, we have 



F R with F 



(26) 








/'2 








1 





1 


1 





-1 








P-2 








(27) 



The identity 

U*RU = F^rfSU (28) 
follows immediately. As U is a reflection in the relative coordinate, we can easily calculate 

/ 1 o o o \ / \ 

f 
-1 
\0 -1 / 



U^SU = GS + B with G 



D 





2a 
\ o / 



(29) 



Using eq. (J25j) and (J2HJ), we arrive at 



U ] RU = F~ 1 GFR + F- l B . 



Thus the covariance matrix a' is related to the covariance matrix a of eq. (j^lj) by 



(30) 



a 1 = F- 1 GF a F T G T 



(31) 



(32) 



We write cr' in the form 

' A C 
C T B 

where A, B and C are 2 by 2 matrices. Then a straightforward computation of the product 
of matrices in eq. (}3*Tj) yields 



A 



C 



2^f + ^ 





1»\ , (A/j) 2 
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2^0-? + 





2^1 (A M ) 2 







where we used the abbreviation 



AfJL = H\ — jl 2 



(33) 



(34) 



Looking at the definitions of covariance matrix and reduced density operator, it is clear 
that the covariance matrix which corresponds the reduced state p red of eq. (j!8|) is given by 



the submatrix A of eq. 



It is known 
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I ~\ ^1 that the von Neumann entropy as 



defined in eq. (J2J) - of a one mode Gaussian state p red is 

S (p rcd ) = (d + log 2 (d+^j-(d- log 2 (d-^j , (35) 



where the quantity d > 1/2 is given as 

(A/i) 2 pj aj (4 °l 



df = det(A) = det(B) = i/iffx* + (A/i) 



(in a 



(36) 



If we insert the (positive) value of d(fjLi, P2, erf, of) into eq. (|33j) . we get the desired asymp- 
totic entanglement in terms of the four parameters pi, p 2 , of, a\. 

Instead of calculating the von Neumann entropy of the reduced density operator one can 
also consider - as done in 1J] - the purity of the state p red 

V{p- d ) :=Tr(y cd ] 2 ) . (37) 

. n 

The smaller V (p J, the larger is the generated entanglement. It is known [9] that the 
purity of a Gaussian state like p red is given by 

V (p rcd ) = jl = -3- . (38) 

Thus the purity is - as the von Neumann entropy - characterized by the determinant of the 
covariance matrix A. 

We now discuss our result for the entanglement as expressed in eqs. fl33j) and (jSHJ). We first 
note that d 2 depends only on the two width ratios o\ja\ 1 (y\ja\ and on the mass fractions 
pi, p2 of the two particles. We see that for a large d - and therfore a large entanglement - 
it is necessary that one of the two width ratios is large. In contrast to that the quantities 
< pi, p,2 < 1 enter in formula (j3T?j) not as a quotient, so that choosing large mass ratios 
does not generate a large d. We will discuss the conditions for the production of a lot of 
entanglement in detail in section HVl 

One can easily check that always d > 1/2 and thus S (p rcd ) > 0. If d = 1/2, then 
the entanglement vanishes according to eq. ()35j) . No entanglement is generated for the two 
conditions (cf. eq. ((23)) 

1. pi = P2 = 1/2, i. e. the two scattered particles have equal mass. This result was 
already obtained in [l^j ]. 

2. pi o\ = P2(y\, i. e. the heavy particle has a small width, the light particle has a large 
width. 
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IV. THE WAVE FUNCTION AFTER SCATTERING 



In this section we study the wave function go (xi, X2) that corresponds to the state \go) 



2l\ . This will give us more insight into the results in eqs. (|35|) and (J36)) for the generated 
entanglement. In complete analogy to the calculation of the co variance matrix a' in the last 
section, we can obtain the wave function of the state |go) = (1 ®Pa) |/o) in three steps from 
the initial wave function fo(xi, x 2 ) of eq. (jHJ): 

1. Replace in fo(xi, x 2 ) the coordinates (x\, x 2 ) with the coordinates (x s , x r ) using the 
inverse to relation 

2. Apply the operator (U ® P a ) by changing the variable x r to (2 a — x r ). 

3. Go back to the coordinates {x\, x 2 ) with relation (JHJ). 
The result of this procedure is 

ga{xu x 2 ) = 4>g (2 \i 2 x 2 + (/ii - \i 2 ) xr, Qi-2[i 2 a, -K, a 2 ) 

x (p G (2 /i! xi - (/ii - /i 2 ) x 2 ; -Q 2 + 2 /ii a, K, of) , (39) 

where the Gaussians <Pg( ) have been defined in eq. Q. Since the entanglement of a 
Gaussian state depends only on its covariance matrix, any displacement of the wave function 
in position or momentum space does not change its entanglement. Thus instead of go(xi, x 2 ) 
in eq. (J39)) we introduce the simpler wave function g(x%, x 2 ) that keeps only the relevant 
quadratic terms in x±, x 2 in the exponent 



K) exp(- 



[(/■*! ~ /fg) Xi + 2fi 2 X 2 } 2 

2a\ 

[2/iiXi - (fii - fi 2 )x 2 } 2 



x a (aj) exp ( - ^ ^ ™"« ) . (40) 
This wave function can be rewritten with a quadratic form in the exponent 



9o(xi, x 2 ) = a (of) a (of) exp ( - ix T L T £Lx) , (41) 

where we have introduced x T := (xi, x 2 ) and the two matrices 

A*-*. 2/i 2 \ ^ = /l/of ( ^ 
2/i! \ 1/of 
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Defining the matrix in the exponent of the wave function in eq. ()41j) as 

M:=L T SL , (43) 

a comparison with eq. ()40|) yields for its entries 

(/ii-/i 2 ) 2 , 4/x 2 . (M1-M2) 2 

M11 ; — 9 1 5- , M22 •— — o- H o , 

a{ cr| erf cr| 

M 12 = M 21 :=2(/ii-/i 2 ) fe-^) • (44) 

^From eq. (|4*Tj) we see that the wave function factorizes into a product of particle 1 and 2 
wave functions, if and only if M i2 = 0. Eq. f)44j) shows us that this is the case when \i\ = fi 2 
or when \i\ a\ = /i 2 of - in accordance with our findings in the last section. 

We now discuss the two cases with factorizing wave functions in more detail. In the case 
Hi = /i 2 the widths of, of of the particles are interchanged after the scattering process, as 
the formulas for Mn and M 22 in eq. f)44j) show 22] . More generally, it can be shown that 
the two particles exchange the shape of their wave functions in the scattering process, if the 
two particle masses are equal. When \ii erf = iii<j\, a short calculation shows that 

M xx = \, M 22 = \, (45) 
a 1 cr 2 

and thus both particles keep their initial width. That no entanglement occurs for fix a\ = 
fi 2 a\ depends crucially on the Gaussian shape of the incoming particle wavefunction. For, 
the two mixed terms proportional to x\ x 2 in the two exponents of eq. (|4()jl have opposite 
sign, if fii a\ = [x 2 cr 2 , and therefore cancel each other. 

Points {x\, x 2 ) where the wave function go(%i, x 2 ) has the same value form an ellipse. 
For example, we may select those points (xi, x 2 ) where the exponent in eq. (j41j) is —1/2. 
This leads to an analytical expression for an ellipse Ep as follows 

Ep : x T Mx=l . (46) 

Generally this ellipse is oblique to the X\, x 2 axes; this is precisely then the case when 
entanglement occurs. The semimajor axis A and the semiminor axis B of the ellipse Ep in 
eq. fJ46jl are determined from the two eigenvalues Ai, A 2 of M 

A = ^rTi B = 4r> (Ai>A 2 ) . (47) 

V A 2 y/ A\ 
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FIG. 1: Here the initial ellipse Ej with semiaxes a±, 02 and the final ellipse Ep with semiaxes 
A, B are depicted. The rotation angle 9 of Ep is approximately 63, 4° and thus the length of the 
semimajor axis A is more than doubled in comparison to a\, whereas the semiminor axis B is 
shrunk by the same factor in comparison to 02. 

Let 9 be the angle between the semimajor axis of Ep and the X\ axis. 9 is determined by 
the eigenvectors of M. In fig.Qwe depict the geometric quantities A, B and 9. We can also 
associate to the initial wave function fo(xi, x 2 ) in eq. (JBJ) an ellipse by keeping - as done for 
the wave function go(xi, x 2 ) - only the quadratic terms in the exponent yielding 



where the diagonal matrix S has been denned in eq. (J42j) . The quadratic form in the 
exponent of equation (J48|) defines an ellipse Ej where the major axes are the widths ai, a 2 
of the two particles. We adopt the convention that a\ > erf. Thus we can understand the 
effect of scattering in a geometric way as a transformation of the initial ellipse Ej into the 




(48) 



12 



final ellipse Ep (see fig. HJ). Because of eq. P3*Jl and detL = —1, 

det M = det E = -3— 2 , (49) 

and thus the of the ellipses Ej and £7f are equal. 

In general we get for the major axes of the final ellipse (cf. eq. (@7j)) rather complicated 
expressions in our parameters /jli, of, o\. A significant simplification takes place, if we 
consider the case that one particle has a much greater width than the other one, i. e. the 
semimajor axis of the initial ellipse Ej is much longer than the semiminor axis (o~i ^> gq). 
This is also the most interesting: only then a lot of entanglement can be created 

via scattering, as we have remarked before. When <j\ 02, it holds approximately that 
the semimajor axis of the ellipse Ej is mapped via the transformation L in eq. ()42|) to the 
semimajor axis of the ellipse Ep. This means that 



P :=L y q j =o, \ - I , (50) 
and the length A of the semimajor axis in Ep is 



A w ||P|| = y/QM o x , (51) 

where we define the quadratic polynomial 

Q(x) := 8x 2 -Ax + 1 . (52) 

Since the areas of the ellipse Ej and Ep are equal, we get for the length of the semiminor 
axis 

B w / 2 . (53) 

In the same approximation (a\ ^> 02) the angle 9 of the semimajor axis (cf. figQ) can be 
read off from the vector P in eq. ()50|) as j3] 

9 = arctan f 2 ^ ) . (54) 

If we analyze these results, we see that for < [i\ < 1 the polynomial Q(fii) in eq. ()52j) has 
values between the minimum 1/2 (for fii = 1/4 ) and 5 (for /ii w 1). Thus according to 
eqs. (I51jl and ()53j1 the semimajor axis A can be elongated by up to a factor a/5 or shortened 
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"0 6 




FIG. 2: Here the quantity d is plotted as a function of \i\ for fixed width ratio c\ja2 = 10. The 
solid line shows the leading term in formula the dahed line the exact value in eq. 

by up to a factor y/2. The semiminor axis is always scaled by the inverse factor. Looking 
at the angle 9 in eq. (JH3|) . we find 



Ml > ^2 
Ml < ^2 



7T 



arctan(2) < 6 < — 

- < 9 < 71 . 

2 



(55) 



If Hi — \ii — 1/2 we have 8 = n/2 which is consistent with the generation of no entanglement. 
If fix — > we get 6 = 7T from eq. (JS1J), and there is no entanglement as well. In our 
approximation the case /ii corresponds to the case /ii a\ = /i 2 o"2 studied above, since 
terms of the magnitude o\ja\ have been neglected. In general, our analysis shows that 
the scattering process leads to more complicated transformations of the ellipse than in the 
situation of squeezed states hitting a beamsplitter where the initial ellipse is only rotated 

(cf. Q). 

How can we understand the circumstances the production of a lot of entanglement in this 
geometric picture? First, for 0\ ^> a 2 we calculate the leading term in ai/a 2 of the quantity 
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d from eq. (|3T)j) 

d ~ |2 /ii — 1| fix — . (56) 

For the fixed value of <J\l<Ji = 10 we plot in fig. |2]the exact d and this approximation as a 
functions of fix and find very good agreement. Thus we will use the simpler form of eq. (JoT)|) 
for our subsequent discussion. 

According to fig. El (sold line) we find a local maximum of d at fix = 1/4 and the absolute 
maximum at fix — 1. The value at the absolute maximum is greater by a factor 8 than that 
at the local one. If one looks only at the angle of the final ellipse Ep, one would expect 
the absolute maximum for fix = 1/4, since for this value we get the angle 9 = 3n/A from 
eq. (j31J). Here the ellipse E F is more tilted in the xx, x 2 plane than at the value fix = 1 
where we obtain 9 = arctan(2) pa 63,4°. The explanation why the value of d for fix — 1 is 
greater than that for /xi = 1/4 lies in the factor a/ Q(fix) in eqs. (|BTj) and (JH3j) that can lead 
to a stretching or shrinking of the initial ellipse Ei . For, at fix = 1/4 the polynomial Q(fix) 
has - as mentioned before - its minimum with value 1/2 and thus 

Ajflx = 1/4) O"! 1 (7! 

5^ = 1/4) = = 1/4) ^ = 2 ^ • (57) 
In contrast to that, for /ii = 1 we get 

f^ = Q( Al = l)^ = 5^ . (58) 
S(/i! = 1) cr 2 cr 2 

Thus the ratio between semimajor and semiminor axis is augmented by a factor 5, whereas 
for fix = 1/4 it is diminished by a factor 2. This effect outweighs the effect of the angle 9 
and explains why the absolute maximum is located at fix — 1. 



V. CONCLUSIONS 



We have studied in detail the scattering of two Gaussian wave packets with hard core 
repulsion. The generated entanglement could be calculated analytically (cf. our central 
result in eqs. and ). The special hard-core potential that we have considered yields 
Gaussians as asymptotic wave functions. Therefore, we could apply results about Gaussian 
states. 

As we vary for the two particles the ratio of their masses and the ratio of their initial 
widths, it turns out that a great width ratio is necessary for the production of much entan- 
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glement in the scattering. If we vary for fixed and large width ratio the mass ratio of the 
two particles, maxima of entanglement production are reached, 

1. If the mass of the particle with large width is much greater than that of the particle 
with small width. This is the absolute maximum. 

2. If the mass of the particle with large width is approximately 1/3 of the mass of the 
particle with small width (local maximum as shown in fig. |2J). 

By contrast, no entanglement is generated, 

1. If the two particles have equal mass. 

2. If for masses and widths of the two particles the relation fi\ o\ = //2 o\ holds, i. e. one 
particle has large width and small mass, the other particle small width and large mass. 

These results about maxima and minima of entanglement are in our opinion not at all 
obvious, even though scattering at a hard core potential is one of the simplest scattering 
processes. 

We have gained additional geometric understanding in our results by looking at the wave 
function after scattering. This Gaussian wave function defines an ellipse whose properties 
determine the degree of entanglement after scattering. This shows strong parallels to studies 
of entanglement in quantum optics where an ellipse can be defined for squeezed light hitting 
a beam splitter. However in contrast to the beamsplitter, the scattering process can not 
only rotate the ellipse but also stretch or squeeze it which is relevant for the amount of 
entanglement. 

Our results indicate that there remains much to be done in the subject of entanglement 
production via scattering. Beyond the special case of hard core repulsion one should look 
at more realistic interaction potentials which allow for reflection and transmission of the 
particles. Also the modification of our results for bosons and fermions poses an interesting 
problem. We leave these tasks to future studies. 
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